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Abstract. Catalyzed big bang nucleosynthesis (CBBN) can lead to an overproduction of 6 Li in 
gravitino dark matter scenarios in which the lighter stau is the lightest Standard Model superpart- 
ner. Based on a treatment using the state-of-the-art result for the catalyzed 6 Li production cross 
section, we update the resulting constraint within the framework of the constrained minimal super- 
symmetric Standard Model (CMSSM). We confront our numerical findings with recently derived 
conservative limits on the gaugino mass parameter mx/2 and the reheating temperature Tr. 

PACS. 12.60.Jv Supersymmetric models - 95.35.+d Dark matter 



1 Introduction 

Big Bang nucleosynthesis (BBN) is a cornerstone of 
modern cosmology that allows us to probe physics 
beyond the Standard Model. It has been realized re- 
cently that the presence of massive long-lived negative- 
ly charged particles X~ at the time of BBN can have 
a substantial impact on the primordial light clement 
abundances via bound-state formation [1-11]. 

In scenarios in which the gravitino G is the lightest 
supersymmetric particle (LSP), a long-lived X~ may 
be realized if the lighter stau t\ is the next-to-lightest 
supersymmetric particle (NLSP). In particular, a t\ 
NLSP can be accommodated naturally in the frame- 
work of the constrained minimal supersymmetric Stan- 
dard Model (CMSSM) [12-14,4, 15] in which the gaug- 
ino masses, the scalar masses, and the trilinear scalar 
couplings are assumed to take on the respective univer- 
sal values mi/2, mo, and A at Mgut — 2 x 10 16 GeV. 
There the stau emerges as the lightest Standard Model 
superpartner in a large part of the CMSSM parame- 
ter space. Since the couplings of the stau to the grav- 
itino are suppressed by the (reduced) Planck scale, 
M P = 2.4 x 10 18 GeV, n will typically be long-lived 
for conserved R-parity 1 and thus r-T = X~ . 

Then t{~ and 4 He can form bound states, ( 4 HeT 1 ~), 
and too much 6 Li can be produced via the catalyzed 
(CBBN) reaction [1] 

( 4 He?r)+D^ 6 Li + ?f . (1) 

This happens at temperatures T ~ 10 keV [1] when 
standard BBN (SBBN) processes are frozen out. The 
observationally inferred bound on primordial 6 Li then 
severely restricts the t\ abundance at those times 2 and 
thereby the t\ lifetime . 
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1 For the case of broken R-parity, see, e.g., [16,17]. 

2 In this work we assume a standard cosmological history 
in which t\ was in thermal equilibrium before decoupling. 



For conserved R-parity, the gravitino LSP is stable 
and a promising dark matter candidate. After infla- 
tion, gravitinos are regenerated [18] in thermal scat- 
tering of particles in the primordial plasma. The re- 
sulting gravitino density X2~ p will contribute substan- 
tially to the dark matter density /2dm if the radiation- 
dominated epoch starts with a high reheating tem- 
perature Tr [19-21]. In addition, gravitinos are pro- 
duced in stau NLSP decays with the respective den- 
sity ^ TP [22, 12, 23]. 3 

In this work we study gravitino dark matter scenar- 
ios within the framework of the CMSSM in which t\ 
is the NLSP. For two exemplary parameter scans, we 
compute = /?~ p + i?~ TP at every point in the as- 
sociated {mo,m 1 / 2 ) planes and compare it with J7dm- 
The Tr x -dependent exclusion boundary on the stau 
abundance from 6 Li overproduction allows us to in- 
fer the cosmologically disfavored CMSSM region. Our 
numerical findings are confronted with the recently de- 
rived conservative limits on and Tr [11]. We dis- 
cuss the present status of the relevant BBN constraints 
at the end of Sec. 4. 



2 Gravitino Dark Matter in the CMSSM 

In the CMSSM, the superparticle mass spectrum is 
determined by specifying mo, mi/2, Aq, the ratio of the 
two MSSM Higgs doublet vacuum expectation values, 
tan P, and the sign of the higgsino mass parameter \i. 
Then either the lightest neutralino Xi or the lighter 
stau T\ with respective masses and m^o is the 
lightest Standard Model superpartner and hence the 
NLSP 4 whose mass is denoted by ttinlsp- 



3 We do not include gravitino production from inflaton 
decays; cf., e.g., [24,25] and references therein. 

4 A stop ti NLSP is not feasible in the CMSSM [26]. 
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With the gravitino LSP, the relic density from NLSP 
decays reads 

Qf p h 2 = m d Y^ SpS (T )h 2 /p c , (2) 

where is the gravitino mass. The quantity ^nlsp = 
n NLSp/ s denotes the NLSP yield where ^nlsp ^ s ^ ne 
number density at decoupling and s — 2ir 2 g«s T 3 /45 
the entropy density; p c /[s(T )h 2 } = 3.6xl0~ 9 GeV [27]. 
We obtain i^LSP by employing the computer program 
micrOMEGAs 3.17 [28] which we feed with the super- 
particle mass spectrum computed with SuSpect 2 . 34. 5 

The upper panels in Fig. 1 show contours of I^lsp 
(solid) and tonlsp (dotted) in the (m 1( / 2 ,mo) plane 
for A = 0, n > 0, (a) tan j3 = 10 and (b) tan (3 = 30. 
Above (below) the dashed line, to^o < m?j (to^ < 
m^o). The medium gray and the light gray regions 
at small 7774/2 are excluded respectively by the mass 
bounds m~± > 94 GeV and TOh > 114.4 GeV from 
chargino and Higgs searches at LEP [27]. For tan/3 = 
30, tachyonic sfermions occur at points in the white 
corner labeled as "tachyonic." 

Let us now explore the parameter space in which 
the relic gravitino density matches the observed dark 
matter density [27] Ol^h 2 = 0.1051^^ where h is 
the Hubble constant in units of 100 km Mpc _1 s _1 . 
Then, Tr and to^ appear in addition to the CMSSM 
parameters. In the lower panels of Fig. 1, the shaded 
(green in the web version) bands are the (7774/2, 777 ) 
regions in which 

0.075 < Qfh 2 + Hf p h 2 < 0.126 (3) 

for the indicated values of Tr and for (c) = too and 
(d) 777^ = 10 GeV. For f2~ p , we use expression (3) of 
Ref. [20] 6 while J7~ TP is computed as explained above. 
In the dark shaded regions at larger too, the gravitino 
is not the LSP. 

3 Catalyzed BBN of 6 Li in the CMSSM 

We now focus on the T\ NLSP region. For typical val- 
ues of Y± c = Vj^£ c SP /2 [see Figs. 1 (a,b)j, the amount 

of 6 Li produced in the CBBN reaction (1) can be as 
high as 6 Li/H|cBBN = 10~ 7 [1,5,11] and hence far in 
excess of the observationally inferred upper limit on 
the primordial 6 Li abundance [30] 

6 Li/H| obs < 2 x 10- 11 . (4) 

F° r T ?i ^ 5 x 10 3 s [1,5,6,8,11], however, the staus de- 
cay before (1) becomes too efficient so that 6 Li/H|cBBN 
can be in agreement with (4). From the -dependent 
upper limit on Y3- c computed in [11] by solving the 



5 We use the following Standard Model parameters: 
m t = 172.5 GeV^ m b (m b ) MS ' = 4.23 GeV, af^(m z ) = 
0.1172, and a^™(ro z ) = 127.90896. 

6 For the definition of Tr, see Sec. 2 in Ref. [15]. 



Boltzmann equations containing the state-of-the-art 
result for the catalyzed 6 Li production cross section [5] , 
we obtain the long-dash-dotted (red in the web ver- 
sion) lines in the lower panels of Fig. 1. The regions 
to the left of theses lines are cosmologically disfavored 
because of overproduction of 6 Li. 7 Note that only the 
constraint from the primordial D abundance on hadronic 
energy release [29,31] in f\ decays [23,13,32] can be 
more severe than the one from catalyzed 6 Li produc- 
tion [4, 33, 15, 7]. 8 This is shown by the short-dash- 
dotted (blue in the web version) lines in panel (c) 
which exclude the region in the T\ NLSP region above 
these lines. In panel (d) the D constraint does not ap- 
pear; for details see [32,15]. Contours of the NLSP 
lifetime are shown by the dotted lines. 



4 Discussion 

The position of the Li constraint is governed by the 
stau lifetime. From the constraint < 5 x 10 3 s, 
conservative limits on the gaugino mass parameter, 

/ 777 ~ \ 2/5 

^/ 2 >0-9TeV( I ^ F ) , (5) 
and the reheating temperature, 

/ 777 ~ \ 1/5 

T ^ 4 - 9xl ° 7GcV (itugv) (6) 

have been derived recently [11]. Note that these limits, 
obtained in the framework of the CMSSM, do only 
depend on mg. 

The limit (5) emerges since scales with TO1/2 
[see Figs. 1 (a,b)] and since is fixed once TOg and 
777^ are specified. The choice = 10 GeV in Fig. 1 (d) 
allows for an immediate comparison of our numerical 
findings with (5) and (6). Only in the vicinity of the 
dashed line, i.e., in the T\-x\ coannihilation region, 
the position of the 6 Li constraint approaches its con- 
servative lower limit (5). This is because T\ becomes 
heavier for larger too which shortens for fixed mg. 
Contrariwise, the splitting between the actual position 
of the 6 Li constraint and (5) is larger for smaller too- 
At Too = 10 GeV, this is more pronounced in Fig. 1 (d) 
than in Fig. 1 (c). This results from the fact that the 
increase in tan (3 leads to a decrease in so that 
becomes larger for fixed TOg. 

The limit (6) relies on thermal gravitino produc- 
tion only, j?~ p ~ Tr [19-21]. Thus the upper limit on 

Tr can only become more stringent by taking i?5 TP 
into account. In Fig. 1 (c) we have fixed to^ = to . 
Thereby, the non-thermal production (2) becomes more 
important for larger values of mo. In addition, Y"^ cc 

takes on its maximum at a given TO1/2 in the Ti-Xi 
coannihilation region. This leads to the bending of the 
bands (3) towards lower TO1/2. From Figs. 1 (c,d) we 

7 In this regard, cf. discussion at the end of Sec. 4. 

8 Additional constraints on the primordial light elements 
from CBBN can be found in [4,6,7,10]. 
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Constraints on the Reheating Temperature in Gravitino Dark Matter Scenarios 



tan = 10, Ao = 0, ji > tan/3 = 30, A = 0, /x > 




(c) m 1/2 [GeV] (d) m 1/2 [GeV] 



Fig. 1. The (m^^mo) planes for Ao — 0, /j, > 0, and the choices (a,c) tan/3 = 10 and (b,d) tan/3 = 30. Above (below) 
the dashed lines, m-o < (m^ < m~o). The medium gray and the light gray regions at small m\/i are excluded 

respectively by chargino and Higgs searches at LEP. In the upper panels, contours of i^xsp (solid) and jtinlsp (dotted) 
are shown. In the associated lower panels, we consider scenarios with (c) mg = mo and (d) mg = 10 GeV. The light, 
medium, and dark shaded (green in the web version) bands indicate the regions in which 0.075 < Q d h 2 < 0.126 for 
Tr = 10 7 , 10 s , and 10 9 GeV, respectively. In the dark gray region, the gravitino is not the LSP. The dotted lines show 
contours of the NLSP lifetime. With the Ti NLSP, the region to the left of the long-dash-dotted (red in the web version) 
line is disfavored by the primordial 6 Li abundance. Constraints from primordial D disfavor the t\ NLSP region above the 
short-dash-dotted lines. The region to the left of the thin gray (pink in the web version) line is disfavored by 3 He/D [29]. 
On the solid vertical line (violet in the web version) mj = 2.5 TeV. 



thus find Tr < 10 7 GeV confirming that (6) indeed 
provides a good conservative estimate. The bound Tr < 
10 7 GeV can be very restrictive for models of inflation 
and baryogenesis. 

Since for a t\ NLSP typically <C m\ , 2 , it is the 
gaugino mass parameter mi/ 2 which sets the scale for 
the low energy superparticle spectrum. Thus, depend- 
ing on m^, the bound (5) implies high values of the 
superparticle masses which can be associated with a 



mass range that will be difficult to probe at the Large 
Hadron Collider (LHC). This is illustrated by the ver- 
tical (violet in the web version) line in Figs. 1 (c,d) 
which shows the gluino mass contour wig = 2.5 TeV. 9 
In the above considerations we have assumed a 
standard thermal history of the Universe during the 
radiation-dominated epoch. A substantial entropy re- 

9 Note that the mass of the lighter stop is mjr — 0.7rrig 
in the considered f\ NLSP regions with mi, > 114.4 GeV. 
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lease after the decoupling of the NLSP but before the 
onset of BBN may dilute Y~iL 0C [34] such that 6 Li/H| CB BN 

respects (4) even for r? , > 5 x 10 3 s [15,5]. The actual 
amount of entropy required, however, is model depen- 
dent. For example, if entropy is released by an out-of- 
cquilibrium decay of a massive particle species <j>, the 
presence of the energy density p<f, during NLSP decou- 
pling can affect 5^£g P . Furthermore, the branching ra- 
tio of the (j) decays into t\ and/or gravitinos may be 
substantial; see, e.g., [24,25]. An illustrative scenario 
taking into account the former effect but neglecting 
the latter can be found in [15]. As argued in [35], how- 
ever, a concrete realization of a large entropy release 
by (j) decays in the narrow time window after NLSP 
decoupling and before BBN might be rather difficult 
to accomplish in the considered scenarios. 

Let us finally comment on the present status of 
BBN constraints on gravitino dark matter scenarios 
with a long-lived charged slepton NLSP. It has recently 
been pointed out in Ref. [10] that that bound-state for- 
mation of X~ with protons at T ~ 1 kcV might well 
reprocess large fractions of the previously synthesized 
6 Li. The drop in the 6 Li abundance at t ~ 3 x 10 6 s 
found in [10] is rather drastic. As can be seen from 
the associated contour in Figs. 1 (c,d), this could 
reopen a 6 Li conform window in the collider friend- 
ly region of low and low mj, provided that the 
LEP Higgs bound is respected. Unfortunately, at the 
time of writing, the existing uncertainties in the rele- 
vant CBBN nuclear reaction rates in [10] make it diffi- 
cult to arrive at a final conclusion on 6 Li/H|cBBN for 
t ti > 3 x 10 6 s. In this region, however, the 3 He/D con- 
straint on electromagnetic energy release [36] becomes 
severe and can exclude > 10 6 s [13,4,7,10]. This 
is illustrated in Figs. 1 (c,d) by the thin gray (pink in 
the web version) line which is obtained from Fig. 42 of 
Ref. [29] with E vis = 0.3 (m\ -m 2 ~ + mj)/2m fl . Be- 
cause of the strong sensitivity of the 3 He/D constraint 
for 10 6 s < r ?1 < 10 7 s [29,31], it still remains dif- 
ficult to decide whether small cosmologically allowed 
islands would exist in the CMSSM parameter space 
for r fl > 5 x 10 3 s. 10 



5 Conclusion 

We have considered gravitino dark matter scenarios in 
which T\ is the NLSP. In exemplary CMSSM scena- 
rios we have demonstrated that our recently obtained 
limit [11] T R < 4.9 x 10 7 GcV(m g /10 GeV) 1 / 5 from 
catalyzed 6 Li production is indeed conservative. In par- 
ticular, taking into account .f2~ TP , the Tr limit can be- 
come considerably more severe. Furthermore, we have 
shown explicitly that the 6 Li constraint can exclude 
nig < 2.5 TeV. The cosmologically favored region can 
thus be associated with a mass range that will be very 
difficult to probe at the LHC. 



With a highly fine-tuned m^-mg degeneracy leading 
to -Evis — > 0, any bound on energy release can be evaded. 
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